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The purpose of the grant entitled "Continuation of the Experimental Investigation of 
Surface Tension Driven Convection as a Feasibility Study for a Microgravity Experiment" 


« 


was the completion of an experiment checking on the properties of supercritical surface 
tension driven Benard convection. Since there is no theory describing the nonlinear aspects 
of surface tension driven Benard convection in a fluid layer of large aspect ratio, or in 
theoretical terms in a fluid layer of infinite horizontal extent, the only way to obtain 
information about this topic was to make an experiment studying this problem. When the 
proposal for this grant was submitted we had already obtained the first results of such 
experiments. They indicated that the size of hexagonal surface tension driven cells in a 
shallow layer of silicone oil decreased with increased Marangoni number or increased 
temperature difference across the fluid layer. That was a startling result since it is well- 
known that the cell size of Rayleigh-Benard convection increases with increased 
supercritical temperature difference. In order for our preliminary results to stand, the 
finding of a decreasing cell size with increased temperature difference had to be made 
certain. Over the time-span of the grant this has been done, the results of our experiments 
are spelled out in detail in the attached manuscript: "The Wave Numbers of Supercritical 
Surface Tension Driven Benard Convection" by E. L. Koschmieder and D. W. Switzer. 
This manuscript has been submitted to the Journal of Fluid Mechanics and has been 
accepted for publication. All detail of the experiments are described in this paper. The end 
result of the experiments supported by this grant is that, indeed, the cell size of supercritical 
surface tension driven hexagonal Benard cells decreases with increased temperature 
difference, and that the cell size is unique , an equally fundamental feature. 


r- 

ro m 

r- (/) s o 

ft TJ <\J 

I — O' 

<M 0-0 

O' CO 

Z D o 


o 


nt 

CD 



a 

O I 

cc 


< 


'Z 

OJ 

CC 

1 

LL 

• 

O 2T 

> 

LU 

• •m 

> 

c 

CL »-» 

z> 

lu or 

CO o 

HI 

X 

CD 

3 Z 

X 

Z C 

d 

•— t 

h- 

LU (/) 


> z 

< LU 

~x t— 

L 

LU LU 

O 

X o 

a 

t- < 

a) 

LL. 

0 c 

CC 

D 

f— 

,0 V) 

fO 

IT 

c 

CO — 1 

•— 

O' < 

u. 

CO O 

rH ►- 

z 

1 H- 

o 

or i-t 

t—t 

u a c 

t— 

l <_) 

u 

< OC 

LU 

to UJ 

> 

< CL. 

z 

Z Z3 

o 

w (/) 

o 


O 

(M 


O 

(/> 

U 


Department of Civil Engineering- Austin, Texas 78712-1076 

Environmental and Water Resources Engineering -(512) 471-5602 -FAX (5 12) 47 1-0592 


A 


Hie Wave Numbers of Supercritical Surface Tension Driven Benard 

Convection 

E. L. Koschmieder and D. W. Switzer 

College of Engineering and I. Prigogine 
Center for Statistical Mechanics 

The University of Texas at Austin 
Austin, TX 78712, U.S.A. 


The cell size or the wave numbers of supercritical hexagonal 
convection cells in primarily surface tension driven convection on a 
uniformly heated plate has been studied experimentally in thermal 
equilibrium in thin layers of silicone oil of large aspect ratio. It has been 
found that the cell size decreases with increased temperature difference in 
the slightly supercritical range, and that the cell size is unique within the 
experimental error. It has also been observed that the cell size reaches a 
minimum and begins to increase at larger temperature differences. This 
reversal of the rate of change of the wave number with temperature 
difference is attributed to influences of buoyancy on the fluid motion. The 
consequences of buoyancy have been tested with three fluid layers of 
different depth. 


1. Introduction 

As is well-known the hexagonal convection cells observed by Benard 
(1900) were actually driven by surface tension gradients as was established 
experimentally by Block (1956) and explained theoretically by Pearson 
(1958). Pearson's linear theory assumes that the convective motions are 
driven exclusively by surface tension forces without a contribution of 
buoyancy, that means that the acceleration of gravity was set zero. Onset of 
surface tension driven convection is determined by the critical Marangoni 
number and the cells are characterized by a unique critical wave number. 
The combined surface tension driven and buoyancy problem was studied by 
Nield (1964). Onset of convection in this case is determined in thin fluid 
layers by the critical Marangoni number and in deep fluid layers by the 



critical Rayleigh number, and the Marangoni number and Rayleigh 
number for onset of convection by surface tension gradients in the presence 
of buoyancy are coupled by a simple equation. The convective motions are 
again characterized by a unique critical wave number. As is well-known 
linear theories do not permit a determination of the preferred planform of 
the convective motions which is hexagonal as all experiments show. 

Nonlinear theories dealing with supercritical surface tension driven 
convection in fluid layer of infinite horizontal extent have been published by 
Scanlon and Segel (1969), Kraska and Sani (1979) and Cloot and Lebon 
(1984). Supercritical surface tension driven convection in small circular 
containers has been studied by Rosenblat, Davis and Homsy (1982). The 
papers dealing with supercritical surface tension driven convection in 
infinite layers arrive at the conclusion that the hexagonal pattern is the 
only stable supercritical pattern and hence preferred. Only Cloot and Lebon 
make a statement about the wave number of the supercritical cells, they 
predict that the wave number of the supercritical motions is larger than the 
critical wave number and is nonunique, in the case that the Rayleigh 
number is zero, i.e. is the absence of buoyancy. On the other hand, the 
wave numbers are predicted to be smaller than the critical wave number 
and nonunique for values of the Rayleigh number 400 < R < 500, that means 
when buoyancy is present but small. 

Little is known from experiments about the variation of the size of the 
hexagonal convection cells in surface tension driven convection when the 
temperature difference or the Marangoni number is increased above 
critical. The best piece of information about this question can be found in 
Benard's (1900) original paper, where it is shown in Fig. 23 that the size of 
the cells in very thin and very wide layers varied with the applied 
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temperature difference from larger cells at large AT, to a minimum size at 
moderate AT, to a somewhat larger size again at the critical temperature 
difference. The concept of the critical temperature difference was not 
known to Benard but is apparent from the fact that the cells disappeared at 
this AT when AT was decreased. The only way Benard could vary the 
temperature difference across the fluid was by permitting the metal block at 
the bottom of the layer to cool down from an initial temperature of 100°C, 
thus AT decreased in his experiments. To our knowledge B6nard's 
observation of the variation of the cell size has not been mentioned in the 
modern literature although his observations were discussed in a review of 
the Benard problem by Koschmieder (1974). Dauz&re (1912) published 
convincing pictures showing that the size of hexagonal convection cells 
increased when the applied temperature difference was increased above the 
maximal AT used by Bernard. This confirms the trend towards larger 
convection cells with increased AT reported by B6nard. Since the cell size is 
inversely proportional to the wave number and since B6nard decreased the 
applied temperature difference, whereas we increase the temperature 
difference, we shall here, in order to avoid confusion later on, state 
Benard' s observation as follows. His observations showed that the wave 
number of the motions increases when after onset of convection the 
temperature difference is increased, the wave number reaches a 
maximum, and decreases substantially when AT is increased further. 

The increasing cell size of surface tension driven convection observed 
by Benard and Dauzere follows the same trend as the increase of the cell 
size in Rayleigh-B6nard convection. This has been observed in many 
experiments which have been reviewed by Koschmieder (1974). The 
increase of the cell size of surface tension driven hexagonal convection cells 
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has recently been confirmed by Cerisier et al. (1987a). The observed 
increase of the hexagonal cell size seems to contradict Benard's observation 
that the cell size had a minimum from which the cell diameters increased 
with decreased AT until they disappeared at the critical condition. We will 
prove in the following that Benard's observation of an initially decreasing 
cell size was correct, and that the decrease of the cell size under slightly 
supercritical conditions is characteristic for the variation of the cell size of 
surface tension driven convection. We will also show that the supercritical 
cell size of surface tension driven convection is unique within the 
experimental error. 

2. The apparatus 

The apparatus is a modified version of the apparatus used by 
Koschmieder and Biggerstaff (1986), henceforth referred to as K & B, for the 
investigation of the onset of surface tension driven convection. The bottom 
of the fluid layer was a chromium coated silicon crystal 12.7 mm thick and 
of 17.8 cm diameter. The surface of the coated crystal was plane to 20 
wavelengths and served as the mirror for the shadowgraph visualization. 
The silicon crystal is an excellent thermal conductor, its thermal 
conductivity is X = 0.37 cal/cm sec °C, that is 1/3 of the thermal conductivity 
of copper ( X = 0.94 cal/cm sec °C). The silicon crystal was heated from below 
by a copper block, 5 cm thick and of 17.8 cm diameter. The copper block was 
heated in turn by a regulated electrical current passing through a 
resistance wire. The arrangement with the silicon crystal and the copper 
block should provide a practically uniform temperature at the surface of the 
crystal, by the excellent thermal conductivities of the materials and the 
thickness of the copper block, in combination with the very poor thermal 
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conductivity of the silicone oil (which is of order of 4-1 0" 4 cal/cm sec °C) on 
top of the crystal. 

The silicone oil layer on top of the mirror was bounded by lucite rings 
of 10.48 cm inner diameter and different thicknesses, ranging from 3 mm 
to 1.2 mm machined with an accuracy of order of 10 -2 mm. The nominal 
thermal conductivity of lucite differs from the thermal conductivity of the 
silicone oil by only about 10%. On top of the fluid layer was a thin air gap of 
0.4 mm depth. The thickness of the air layer was held as low as possible in 
order to suppress motions of the air on top of the fluid. The air was 
contained above by a sapphire lid of 5 mm thickness and 13.34 cm diameter 
whose top and bottom surfaces were polished and optically plane. This 
sapphire was used before in the experiment of Koschmieder and Pallas 
(1974), in which such a sapphire was used for the first time in a convection 
experiment. The purpose of the sapphire is to provide a transparent lid of 
good thermal conductivity, which is A, = 0.088 cal/cm sec °C according to the 
manufacturer (Union Carbide). The thermal conductivity of the sapphire is 
about 200 times better than the thermal conductivity of the silicone oil. The 
sapphire was cooled from above by water of constant temperature (± 10' 2 °C) 
circulated at a rate of about 60 cc/sec. By the uniform temperature of the lid 
and by the small depth of the air gap the temperature of the air on top of the 
fluid was made practically uniform, as is required by theory. The sapphire 
was held in place by a lucite ring surrounding the ring containing the 
fluid. The essential features of the apparatus described above are shown in 
Fig. 1. 

The experimental setup complies in a very good approximation with 
the ideal assumptions made usually in theory. These are uniform 
temperatures on a plane perfectly conducting bottom as well as uniform 
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temperatures on top of the fluid layer, which is, however, in surface tension 
driven convection necessarily of poor thermal conductivity because air is a 
very poor thermal conductor. Another assumption made usually in theory 
is that the fluid layer is of infinite horizontal extent, which assumption 
cannot be fulfilled in an experiment. This condition has been approximated 
_by the large aspect ratios of the fluid layers, the maximal value of which 
was 87. The aspect ratio is defined as the ratio of the diameter of the fluid 
layer divided by the fluid depth. The consequences of the presence of the 
lateral boundary were always noticeable. 

The temperature differences applied to the fluid layers are expressed 
in nondimensional form by the Marangoni numbers 

/dS\ ATd 

M ~ (dTj pvic (1) 


and the Rayleigh numbers 


R _ ag_ATd3 (2 

V K 

where dS/dT is the variation of the surface tension coefficient S with 
temperature, AT the temperature difference applied to the fluid layer, d the 
depth of the fluid, p the density of the fluid, v the kinematic viscosity, k the 
thermal diffusivity and a the expansion coefficient of the fluid, with g being 
the acceleration of gravity. The properties of the two fluids used are listed 
in Table 1. The onset of convection of purely surface tension driven 
convection under an insulating surface takes place at Me = 79.61 according 
to Nield (1964) with the critical wave number being a c = 1.993; the onset of 
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purely buoyancy driven convection under an insulating surface takes place 
at R c = 669.0 with the critical wave number being a c = 2.086, according to 
Nield too. We note that according to Pearson, as well as Nield, the critical 
wave number changes when the upper surface is not insulating, and can 
increase by 50% if the upper surface is conducting in the surface tension 
driven case. In the experiments to be described the upper surface condition 
is nearly insulating, that means that the Biot number is near zero. 

The actual temperature differences applied to the fluid layer has to be 
calculated from the temperature difference between the copper bottom plate 
and the temperature of the cooling water on top of the sapphire because no 
temperature sensor can be placed in the narrow air gap or on the fluid 
surface itself without introducing a disturbance. The temperature 
difference across the fluid follows from the formula 

AT 

AT fl = — (3) 

A fl AZ a 
1 + — 

^ A a AZfi^ 

where AT is the temperature difference between the top of the sapphire and 
the copper block, A fl the mean thermal conductivity of the fluid averaged 
over the temperature difference applied to the fluid, A a the mean thermal 
conductivity of the air and AZfl and AZa are the depth of the fluid layer and 
the depth of the air. After the onset of convection the vertical temperature 
distribution in the fluid is no longer linear, that should, however, have no 
influence on the value of A . Of great importance for the determination of 
ATfl is the contribution of the air layer, because the thermal conductivity of 
air (A = 6.1-1CT 5 cal/cm sec °C at 20°C) is so poor. More than 1/2 of the 
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temperature decrease between the sapphire and the copper plate occurs in 
the air layer, although it is so thin. The contribution of the silicon crystal 
and the sapphire to the temperature decrease are negligible. The 
temperature difference ATf] determined this way is not known with an 
accuracy better than ± 5%, mainly because the depth of the fluid layer is not 
known better than 1%, and the depth of the air layer is not known better 
than 2%. 

3. The experiments 

The first series of experiments were made with a fluid layer of 1.90 
mm depth. It became immediately apparent that very regular pattern 
formed, similar to those in Fig. 2, and that the number of cells increased if 
the temperature difference was increased in the slightly supercritical 
range. An increase of the number of cells is equivalent to a decrease of the 
cell size or an increase of the wave number of the motion. Since in 
Rayleigh-B6nard convection the wave number of the convective motions 
tends to be larger than the critical wave number if the heating is time- 
dependent it was decided to eliminate the possible consequences of time 
dependence by heating extremely slowly. This is also of importance for the 
pattern formation because the fluid needs a long time to establish 
equilibrium in the horizontal plane. This time is characterized by the 
horizontal relaxation time D 2 /k, which is of order of 28 hrs. In all 
experiments described in the following the heating was increased in two 
steps during a day so that the temperature increase across the fluid was 
about 0.5°C per day. That meant that the experiments had to be rim 
continuously for a couple of weeks in order to reach a maximal 
temperature difference of about 30°C across the fluid and the air gap, which 
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temperature difference was deemed to be safe for the sapphire and the 
silicon crystal. 

The series of quasi-steady experiments was first made with a lateral 
ring of 1.90 ± 0.01 mm depth, filled to the rim with silicone oil of 1.00 
cm 2 /sec viscosity. Filling the layer to the rim of the ring eliminates 
undesirable effects of a meniscus. The depth of the fluid was, however, not 
known to an accuracy better than ±0.02 mm, because of the accuracy with 
which the fluid depth can be measured with a micrometer. Figs. 2 and 3 
show shadowgraph pictures obtained with such a layer, and Fig. 4 shows 
the increase of the number of cells of the layer when the temperature 
difference was varied. Each point in Fig. 4 originates from a shadowgraph 
picture taken 12 hours after the last increase in heating. The number of 
cells is proportional to the wave number a of the motions, which for 
hexagonal cells is given by a = 4n:d/3L, where L is the side length of the cells 
and d is the depth of the fluid. The wave number as given by the number N 
of the hexagonal cells and the area A of the fluid layer follows from a = 

4rcdV 2.598 IN / 3 VaT According to this formula the wave number of the 
fluid motions is proportional to the number of cells in the fluid layer, and 
the size of the cells A/N is inversely proportional to the wave number. 

The onset of convection was observed at AT C = 7.1±0.1°C. The onset is, 
however, not an unambiguous spontaneous event, the shadowgraph 
picture is initially weak with a substantial background of diffuse white 
light. It is not possible to determine the critical condition with an accuracy 
better than ± 0.1°C. The white lines outlining the hexagonal cells in Fig. 2 
indicate cold fluid, because cold fluid focuses the parallel light with which 
the mirror bottom of the fluid is illuminated from above. The shadowgraph 
pictures have been discussed in detail by Jenkins (1988), the hexagonal 
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pattern in Fig. 2 agrees with Fig. 5a calculated by Jenkins for hexagonal 
cells. The shadowgraph technique is the most suitable visualization for our 
experiments because it permits to do the experiments with the clear fluid 
whose surface tension coefficient is not altered in any way. Noteworthy in 
Fig. 2 are also the consequences of the lateral wall on the cells near the rim. 
These effects can, of course, not be eliminated entirely. The fact that the 
cells at the rim apparently extend beyond the rim seems to be an optical 
illusion brought about by refraction of the light in the fluid adjacent to the 
wall, where the uniformity of the temperature is altered by the different 
thermal conductivities of the silicone oil and the lucite ring. This effect 
increases as AT is increased. 

Increasing AT does not only increase the number of cells and the 
contrast in the shadowgraph, but also changes the shadowgraph picture of 
the hexagonal cells, see Fig. 3. The first sign of change was the appearance 
of bright dots at the vertices of the cell boundaries, as in Fig. 3. These dots 
gradually increase in size when AT is increased and the lines connecting 
the vertices increase gradually in width, although the dots at the vertices 
remain individual dots of greater brightness. We believe that these changes 
reflect the changes of the temperature field in the cells as the conditions 
become increasingly nonlinear. 

When AT approached the maximal AT used in our experiments the 
number of cells seemed to become independent of AT, see Fig. 4. When the 
maximal AT applied was reached the temperature difference was held 
constant for two days, and afterwards AT was decreased just as slowly as it 
had been increased, in order to see whether there is hysteresis in the wave 
number of the motions. In general the number of cells in the layer then 
decreased as is shown in Fig. 4, but the points for decreasing AT do not 
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coincide with the points for increasing AT. There is usually a small gap 
between corresponding cell numbers. Right after the decrease of AT was 
begun the number of cells actually increased with decreased AT. These 
fluctuations of the number of the cells seem to be caused by wall effects. 

The cells which come into the fluid layer when AT is increased do not form 
as individual additional cells in the interior of the layer, but form in small 
groups at the rim, where the largest nonuniformity of the vertical 
temperature gradient is. Likewise the cells that disappear from the layer 
when AT is decreased do not tend to disappear individually. The interaction 
of the cells with the rim is different, depending on whether AT is increased 
or decreased. The difference in the number of cells for increasing or 
decreasing AT at a given temperature difference is therefore not a 
consequence of hysteresis of the wave number of the motions, but the result 
of the boundary conditions at the rim. This will become even more obvious 
when we discuss the 1.2 mm deep layer. 

In order to check on the reproducibility of the cell number as a function 
of AT the same fluid layer whose wave number changes are plotted in Rig. 4 
was heated up again after the pattern had disappeared at AT C and had been 
held at a temperature difference three degrees smaller than AT C for two 
days. The pattern reappeared at a slightly smaller AT C with a slightly 
smaller number of cells as shown in Fig. 4, but these differences result 
from the difficulties pinpointing exactly the critical condition through the 
shadowgraph. Increasing AT again in the same slow way used before 
produced a much more unsteady increase of the number of cells as 
compared to the results of the first temperature increase. There was now 
in particular an extended range in AT over which the cell number 
remained constant. But in the end very near the maximal AT the number 
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of cells was practically the same as it was the first time. The differences of 
the number of cells observed at the same AT after the first and second 
temperature increase only indicate the experimental uncertainty in the 
determination of the number of the cells. When the temperature difference 
was increased further after the uppermost point of the second curve was 
taken the fluid made contact with the sapphire lid. Virtually all 
experiments ended this way. When this happens the fluid apparently 
creeps across the lateral wall of the layer, up the ring carrying the sapphire 
and much of the fluid is then pulled by adhesion to the lid. 

The qualitative results of this experiment are obviously that the 
number of hexagonal convection cells in this fluid layer increased as the 
temperature difference was quasi-steadily increased in the slightly 
supercritical range, that means that the wave number of the convective 
motions increased. It appears that the increase of the number of cells 
comes to an end at the upper end of the curves. It also appears that the 
number of cells or the wave number of the motions is independent of the 
initial condition, within the error of the experiment. These results are in 
agreement with the outcome of two other experiments increasing and 
decreasing AT quasi-steadily with the same fluid in 1.90 mm deep layers. 

The onset of convection occurred in this fluid layer at a temperature 
difference AT = 7.1°C across the fluid. This corresponds to a Marangoni 
number M = 61.8, using the material properties listed in Table 1. The 
Rayleigh number at the onset of convection was 40. The Marangoni 
number Mtheor for onset of convection of the combined surface tension - 
buoyancy problem follows from the relation 
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M_ R_ 

M c + R c = 1 


( 4 ) 


according to Nield (1964). It follows that in this layer Mtheor = 74.8, and 
therefore it is M ce xp = 0.83 Mtheor- The implication of this will be discussed 
later. The wave number of the motions in the hexagonal cells follows from 
the area of the fluid layer, the depth of the fluid and the number of the cells, 
using the relation a = 4jtd/3L, as stated before. With 185 cells at onset of 
convection we find ac = 1.85, while the theoretical critical wave number 
under an insulating surface is ac = 1.993, according to Nield (1964). The 
wave number found will be compared with the wave numbers found in the 
other layers and with the theoretical critical wave number later on. 

The qualitative results with the 1.9 mm deep layer raised the question 
whether they do not contradict the observations of Dauzbre (1912) and 
Cerisier et al. (1987a). On the other hand they confirm the first part of 
Benard's (1900) curve for the number of cells versus temperature 
difference. Since the increase of the number of cells in our experiment 
seemed to have come to a halt at the highest AT we worked with, it was 
decided to extend these experiments to higher Marangoni numbers. This 
can be accomplished by using a fluid of smaller viscosity. The results of two 
quasi-steady experiments made with silicone oil of 0.50 cm 2 /sec viscosity 
and 1.90 mm depth are shown in Fig. 5. It is apparent from this figure that 
the initial increase of the number of cells changes to a gradual decrease of 
the number of cells if the Marangoni number is made sufficiently large. 
This then is in agreement with the observations of Dauzere and also of 
Cerisier et al. in the case of large temperature differences or Marangoni 
numbers. Furthermore Fig. 5 is in complete agreement with Benard's 
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observations, it shows an initial decrease of the cell size, a minimum of the 
cell size, and a subsequent increase of the cell size. 

These observations raise the question how it can come about that the 
variation of the wave number with temperature difference changes sign at 
a moderately supercritical temperature difference. There is, at present, no 
theory that would predict or explain this. We note however that the wave 
number of buoyancy driven Rayleigh-Benard convection decreases with 
increased Rayleigh number. We therefore assumed that the decrease of the 
wave number of the hexagonal cells in the upper part of Fig. 5 is induced by 
buoyancy, which is also present in our experiments because the Rayleigh 
number is not zero, and the value of the Rayleigh number increases with 
increased AT. A way to test this hypothesis experimentally is to vary the 
relative importance of the Rayleigh number, which can be easily done by 
varying the depth of the fluid layer. Since the Rayleigh number is 
proportional to d 3 while the Marangoni number is proportional to d, the 
importance of the Rayleigh number increases in deeper fluid layers. If the 
increase of the wave number of the hexagonal cells is caused, as assumed, 
by the contribution of buoyancy to our experiments, then the decrease of the 
wave number should appear sooner in deeper fluid layers and later, that 
means at higher Marangoni number, in thinner fluid layers. 

The importance of the surface tension effects should increase if the 
fluid depth is decreased. We have therefore determined the variation of the 
number of hexagonal cells in a 1.20 mm deep layer of 10.48 cm diameter. 
Since the critical temperature difference necessarily increases as the fluid 
depth is decreased we had to use for this experiment silicone oil of 0.50 
cm 2 /sec viscosity. The aspect ratio for this layer was 87. Excellent pattern 
with 450 and more cells appeared, as shown in Fig. 6. When AT was 
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increased, the number of cells increased, as in the 1.9 mm deep layer. The 
variation of the number of the cells with increased AT is shown in Fig. 7. 

In this experiment the temperature difference was quasi-steadily 
increased, later on quasi-steadily decreased to subcritical values, and then 
quasi-steadily increased again. In this layer with its increased aspect ratio 
the consequences of the lateral wall should be smaller than with the 1.9 
mm deep layer. The influence of the aspect ratio on the critical wave 
number of hexagonal cells has been investigated experimentally by Cerisier 
and Zouine (1989) who found that within the experimental error the critical 
wave number is independent of the aspect ratio. They also found the 
critical wave number at around 2.1, whereas the theoretical value is 1.993. 
The pattern in the 1.2 mm deep layer was noticeably affected by the lateral 
ring, forming long lines of cells along the ring, instead of individual cells 
adjusting their form to the shape of the wall. When additional cells came 
in as AT was increased they did not come in individually in the interior of 
the layer, but in groups or even lines at the lateral wall. The addition of 
entire lines of cells to the pattern was particularly noticeable when the layer 
was heated for the second time. This was reminiscent of the cell changes 
with the second heating in Fig. 4. Overall the scatter of the data in Fig. 7 is 
about ± 3% and is smaller than the scatter of the data in Fig. 4. 

As can be seen in Fig. 7 the number of the cells increased steadily over 
the range of AT we could work with. The data points end at the same 
maximal temperature difference AT = 30°C across the fluid layer and the 
air gap at which the experiments in Fig. 4 ended. Within the error of the 
experiment the number of cells appeared to be independent of the initial 
condition. This experiment shows that the increase of the number of cells 
(or the increase of the wave number) is more pronounced when the 


15 



contribution of buoyancy to the convective motions is reduced in a less deep 
layer, or when the Rayleigh number is much smaller. The results of this 
experiment confirm the qualitative conclusion reached earlier that the cell 
size of surface tension driven hexagonal convection cells decreases in the 
slightly supercritical range, as long as there is no contribution of buoyancy 
to the motions. 

Onset of convection in the 1.20 mm deep layer occurred at a 
temperature difference of 5.48°C at a Marangoni number M = 61.4 and a 
Rayleigh number R = 18.4. The experimentally determined Marangoni 
number is therefore M ceX p = 0.79 Mtheor- The wave number of the pattern at 
critical with 450 cells is a ce = 1.85, the theoretical critical wave number for 
an insulating upper surface is a c = 1.993. The slope of the variation of the 
wave number with the Marangoni number is Aa/AM = 0.00425, assuming 
that the variation is linear over the entire range observed. We note that the 
variation of the wave number observed in our experiments does not depend 
on the Marangoni number only, but also on the Rayleigh number, because 
as the Marangoni number increases the Rayleigh number increases as 
well, because both numbers depend on the same AT. Since the wave 
number a = a(M,R) it follows that the da/dM we observe is given by da/dM = 
9a/9M + dR/dM-9a/9R. 

The variation of the wave number of hexagonal convection cells was 
also determined in a 3.00 mm deep layer of silicone oil of 1.00 cm 2 /sec 
viscosity of 10.16 cm diameter. Increasing the fluid depth decreases the 
aspect ratio, in this case to nearly 30. That reduced of course the number of 
cells significantly. The pattern remained however clearly hexagonal as 
shown in Fig. 8. The variation of the number of cells with temperature 
difference is shown in Fig. 9. Except for a very small decrease of the 
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number of cells near the onset of convection the number of cells remained 
virtually constant as AT was increased. The initial decrease of the number 
of cells was confirmed by another experiment, which ended when the fluid 
made contact with the lid. In an earlier experiment Koschmieder (1967) 
found that the number of hexagonal cells in a fluid layer of 4.3 mm depth 
and 20 cm diameter increased from 123 to 126 when AT was increased to 
about 2.4 times the critical temperature difference. In this experiment the 
Rayleigh number at onset of convection was about R=230, that means that 
buoyancy probably contributed to instability. The experiment with the 3 mm 
deep layer shows that at this depth the initial increase of the number of the 
cells had nearly disappeared, the experiment with the 4.3 mm deep layer 
shows that the fluid does not pick up the decrease of the wave number of 
buoyancy driven convection very quickly when the fluid depth is increased. 

Onset of convection in the 3 mm deep layer occurred at a temperature 
difference of 5.10°C which corresponds to a Marangoni number M = 71. 

The Rayleigh number at onset of convection was 116, the relation between 
the experimental Marangoni number and the theoretically predicted 
Marangoni number Mtheor is then M cexp = 1.08 Mtheor- The wave number of 
the pattern consisting of 68 cells is a c = 1.86. The critical wave number 
according to theory is a c = 1.99 if the upper surface is insulating. 

We summarize the results of experiments with the different fluid 
layers in Table 2. As we see there the observed critical wave numbers of the 
experiments are in the range from 1.85 to 1.90. The experimental 
uncertainty of the wave numbers is 2-3%. It seems to be fortuitous that the 
value of a cex p is about 1.85 in three of the four layers. We attribute the 7% 
difference between the experimental value 1.85 and the theoretical value 
1.99 to the effects of the lateral wall. A smaller wave number means larger 
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cells, and the cells at the wall are obviously larger than cells in the interior 
of the layer. We also note that a very large fraction of the number of cells is 
affected by the wall; in the 1.90 mm deep layer more than 20% of the cells 
are in contact with the wall. We therefore believe that the observed wave 
numbers at the onset of convection are in reasonable agreement with the 
theoretical prediction. 

The results concerning the onset of convection in the different fluid 
layers are presented in Table 2 in the column with the ratio M ce /M c t- There 
seems to be a decrease of this ratio as the depth of the fluid is decreased. 

But we note that we cannot expect that the experimental Marangoni 
numbers are more accurate than ± 10% on an absolute scale. This is so 
because we do not know the value of dS/dT with an accuracy better than ± 
5%, and we also do not know the value of AT with an accuracy better than a 
couple of percent. This is so because AT has to be calculated with equation 
(3) in which AZ^ is only known to a couple of percent, because the fluid 
depth is only known to about a thousandth of an inch (= 0.025 mm), which is 
a large fraction of the nominal depth of the air gap which is 0.4 mm. The 
experimental Rayleigh numbers are for similar reasons of equal 
uncertainty. On a relative basis, however, the ratios M ce /M c t should be 
much more accurate, because dS/dT drops out, so the observed decrease of 
M ce /M c t seems to be fairly certain. 

However, theory states that there should be no motion in the fluid at 
Marangoni numbers smaller than M c t. But we have observed onset of 
surface tension driven convection at values below M c t before in K & B. Our 
current observations confirm to a surprising degree the experimental 
Marangoni numbers M 2 at which onset of convection in the form of 
hexagonal cells was found in K & B. The onset of convection in a 2.57 mm 
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deep layer of 1.00 cm 2 /sec silicone oil occurred then at M 2 = 76.8, in a 1.81 
mm deep layer of 1 cm 2 /sec silicone M 2 was 61.0, and in a 1.34 mm deep 
layer of 0.50 cm 2 /sec silicone oil it was M 2 = 68.0. Although the depths of the 
layers are not the same as those here, the Marangoni numbers M 2 are 
similar for similar depth in Table 2. There is certainly the same tendency 
for the ratios M ce /M c t to decrease with decreased depth. There is also the 
same apparent failure to continue the decrease of M ce between 3 mm depth 
and 1.90 mm depth to the 1.20 mm layer that we see in Table 2. However, 
the viscosity of the fluid in the 1.20 mm layer was only 0.50 cm 2 /sec, and so 
was the viscosity on the 1.34 mm deep layer in K & B. We do not need to 
expect complete agreement with the results of K & B, because the heating in 
our present experiments was much slower than in K & B and the 
visualization was different, and there are some subtle differences in the 
apparatus. 

Finally we note that in our present experiments onset of convection 
occurred at M ce = 1.08 M c t in the 3 mm deep layer. According to theory 
onset of convection must occur at M c t- The difference between experiment 
and theory may simply be due to the experimental uncertainty in the 
determination of the absolute value of M. On the other hand, the difference 
in the critical value of M may be genuine, and may be caused by the upper 
boundary condition, which is not necessarily completely insulating, as is 
the case when we set M c t = 79.6. The thin air layer on top of the fluid which 
is covered by the well conducting sapphire is not a good representation of an 
insulating upper surface. We have determined the value of the 
nondimensional parameter L which describes conduction on top of the fluid 
layer in K & B, which applies also for our present experiments because the 
setup is essentially the same. We found that L is of order of 0.2 under these 
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conditions. With L = 0.2 the critical Marangoni number is M c = 87 
according to Nield. Our experiments with the 3 mm deep layer are 
compatible with theory if we use M c = 87 for the onset of convection. 

Finally we want to discuss a characteristic defect of the cellular 
patterns which appeared in our experiments. Over the course of the 
experiments the patterns were usually very regular, with truly hexagonal 
cells dominating, except of course right at the wall. Naturally 
imperfections of the patterns occurred. Such imperfections are referred to 
as defects. Defects have to be characteristic in order to be significant. A 
photograph with a defect characteristic for hexagonal cell patterns is 
shown in Fig. 10, where one can see what we call a "star" in the lower part 
of the layer. Such a star can also be found on Fig. 3c of Cerisier et al. 
(1987b), where many other defects are discussed. A star is the combination 
of six deformed cells which have a common center point, each of the six 
cells having only five vertices. There seems to be a point-like disturbance at 
the center of the star, producing a sink. Usually the stars appeared 
spontaneously in the layer after a temperature increase, and usually the 
stars stayed at a fixed position for quite a while, a couple of days or a couple 
of temperature increases. In one particular experiment there were two 
stars simultaneously in the layer. We have noticed the existence of such 
stars much earlier, see Fig. 4 in Koschmieder (1974). In the present series 
of experiments we cannot blame a scratch on the bottom for the formation of 
the star, as we did then. It can, of course, not be excluded that there was an 
impurity at the center of the star, such as a grain of dust or a piece of lint, 
but we have never been able to locate such a disturbance after the 
experiment. It simply may be that the formation of the star is a natural 
way of the pattern to accomplish a change of the wave number. 
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4. Conclusions 

The principal result obtained from these experiments is the 
observation that the wave number of surface tension driven hexagonal 
convection cells increases when the Marangoni number is increased in the 
slightly supercritical range. This observation confirms the same 
observation made by Bdnard (1900) in his original experiments. The 
decrease of the cell size of the hexagonal cells has implications for 
nonlinear stability in general. The increase of the wave number of surface 
tension driven- convection is in startling contrast to the experimentally well 
established decrease of the wave number of buoyancy driven Rayleigh- 
Bdnard convection under supercritical conditions which was discovered by 
Koschmieder (1966). The decrease of the wave number of supercritical 
Rayleigh-Benard convection is discussed in detail in Koschmieder (1974). 
The variation of the wave number of the convective motions is a 
consequence of nonlinear effects and is of prime interest for the nonlinear 
theory of these instabilities. 

The second result obtained from our observation is the apparent 
uniqueness of supercritical surface tension driven convection. This is 
important because according to linear theory the flow can be nonunique at a 
given supercritical Marangoni number, the theoretical nonuniqueness 
permitting wave numbers smaller as well as larger than the critical wave 
number. The nonlinear theory of Cloot and Lebon (1984) also has 
nonunique solutions. We cannot compare our finding of apparent 
uniqueness of the wave number of the hexagonal cells with corresponding 
results of buoyancy driven Rayleigh-Bdnard convection, because it has not 
been established whether the wave number of supercritical Rayleigh- 
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Benard convection is unique of nonunique. However, there seems to be a 
consensus that time-independent supercritical convective motions have 
always wave numbers smaller than the critical wave numbers, and that 
not all wave numbers smaller than a c permitted by the Eckhaus (1965) 
diagram are actually stable. Whether or not the remaining interval of 
supercritical convective wave numbers smaller than a c is only a line plus- 
minus the experimental uncertainty has so far eluded experimental work. 

The third result of our experiments was the observation that the wave 
number of the hexagonal convection cells decreases when the Marangoni 
number exceeds a certain moderately supercritical value. This is in 
agreement with Benard' s findings, as well as those of Dauzere (1912) and 
Cerisier et al. (1987a). There is the possibility that the wave number of the 
cells cannot continue to increase at its original pace because it would 
intersect the marginal curve of linear theory for surface tension driven 
convection. But that does not necessitate that the wave number begins to 
decrease with increased AT. We assume that the decrease of the wave 
number at larger Marangoni numbers is caused by contributions of 
buoyancy to the flow. We find that assumption confirmed by our 
experiments with fluid layers of different depth. Whether buoyancy does 
indeed cause the ultimate decrease of the wave number of the hexagonal 
cells remains to be clarified by theoretical studies or by a microgravity 
convection experiment. 

We thank Mr. C. A. Bielek and Mr. M. Shamieh for their help with the 
experiments. Financial support of this work through the National 
Aeronautics and Space Agency is gratefully acknowledged. A preliminary 
report on part of this work has been given at the IUTAM Symposium in 
Nice 1990 and can be found in Koschmieder (1991). 
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Figure Captions 


Fig. 1. Schematic drawing of the essential part of the apparatus. 

Fig. 2. Shadowgraph picture of the cellular pattern in the 1.90 mm 

deep layer of 1.00 cm 2 /sec viscosity at 1.24 AT cexp . 200 cells. 

Fig. 3. Shadowgraph picture of the cellular pattern in the 1.90 mm 
deep layer of 1.00 cm 2 /sec viscosity at 1.50 AT cex p- 217 cells. 

Fig. 4. The variation of the number of cells in the 1.90 mm deep layer 

of 1 cm 2 /sec viscosity. Open circles are from the first increase 
of AT, full circles are from the following decrease of AT, 
crossed out circles mark the second increase of AT. 

Fig. 5. The variation of the number of cells in the 1.90 mm deep layer 
of 0.50 cm 2 /sec viscosity. Open circles mark the first increase 
— of AT, full circles mark a second experiment with increasing 

AT. 

Fig. 6. The cellular pattern in the 1.20 mm deep layer at 1.32 AT cex p- 
Note the lines of cells along the curved lateral wall. 511 cells. 

Fig. 7. The variation of the number of cells in the 1.20 mm deep layer. 

Open circles mark first temperature increase, full circles the 
following decrease of the temperature difference, and crossed 
out circles mark the second temperature increase. 

Fig. 8. The cellular pattern in the 3.00 mm deep layer at 1.52 AT cexp . 

Fig. 9. The variation of the number of cells in the 3.00 mm deep layer. 

Open circles from an experiment with increasing AT, full 
circles from another experiment increasing AT, which ended 
when the fluid made contact with the lid. 

Fig. 10. A defect on the cellular pattern in the 1.20 mm deep layer in 

the form of a star in the lower center part of the picture at AT = 
1.11 AT cexp - 476 cells. 
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Table 1 


V 

P 

a 

X 

K 

S 

dS/dT 

cm 2 /sec 

g/cm 3 

oc-i 

cal/cm sec°C 

cm 2 /sec 

dyne/cm 

dyne/cm 0 

1.00 

0.968 

0.00096 

3.7-10- 4 

0.001095 

19.36 

-0.050 

0.50 

0.960 

0.00104 

3.6-10- 4 

0.001025 

18.35 

-0.047 


Properties of the fluids at 25°C 


Table 2 


Depth 

mm 

Viscosity 

cm 2 /sec 

a ce 

M C e 

R 

Mce/Mct 

Aspect 

ratio 

3.00 

1.00 

1.855 

71.0 

115 

1.08 

34 

1.90 

1.00 

1.854 

61.8 

40 

0.83 

55 

1.90 

0.50 

1.904 

72.1 

54 

0.99 

55 

1.20 

0.50 

1.851 

61.4 

18 

0.79 

87 


Results of the experiments 
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